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In this talk I report on recent developments and results relevant for LHC phenomenology at next- 
to-leading order QCD. Feynman diagrammatic and unitarity based methods have both seen con- 
siderable improvements and new ideas recently. Current approaches point towards automated 
evaluation of one-loop multi-particle amplitudes. Many results for notoriously difficult processes 
are under construction by several groups and, given the enormous recent progress, it can be con- 
cluded that LHC phenomenology at full next-to-leading order level will become the standard 
approximation soon. 
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LHC phenomenology at NLO QCD 



1. Introduction 

With the start of the Large Hadron Collider (LHC) at CERN we enter exiting times in particle 
physics. Once the machine is fully operating we will be able to explore an energy range which has 
never been reached in a laboratory. Many people in this audience are involved in the enormous 
challenge to master the huge amount of data which will be produced during collisions. I will focus 
in this talk on the theoretical challenge of predicting LHC phenomena as precisely as possible using 
our good old Standard Model. We all very much hope for a discrepancy between well-understood 
data sets and sufficiently precise theoretical predictions. Such a discrepancy will be a signal of 
new and maybe completely unexpected phenomena which will help us to push our understanding 
of nature. 

What should be expected from the LHC? The first and most important issue is the test of the 
Higgs mechanism within the Standard Model which predicts the prominent Higgs boson. The LEP 
experiments set a lower bound of 1 14.4 GeV on its mass from direct searches. As a Higgs boson 
with a mass higher than about 200 GeV would induce radiative corrections which are in conflict 
with electroweak precision measurements, the Standard Model Higgs boson is pretty much nailed 
down to a low mass window. Ironically, this window is not as easy to close as the mass range 
above the two Z boson threshold of about 1 80 GeV, where two muon pairs would provide a gold- 
plated discovery mode for the Higgs boson. The identification of any signal with the Standard 
Model Higgs boson will also entail detailed measurements of quantum numbers like spin and CP 
properties [jl|, ^]. This necessary homework might take even longer than the time we need to find 
new signals (if we are lucky, of course!). The Standard Model has some deficits, e.g. it does not 
contain a natural dark matter candidate, it has no direct relation to gravity and it is sensitive to 
quadratic renormalisation effects. Although the latter is in principle not a fundamental problem 
for a renormalisable theory like the Standard Model, a popular extension is supersymmetry (susy) 
which would not only cure the problem above a sensible susy breaking scale of about 1 TeV, but 
would even allow to view the SM as a low energy limit of some string theory 1 . Supersymmetric 
extensions of the Standard Model have been extensively studied in the last years, for an overview 
see [EL [I], Bt]. Other extensions embed the SM in higher dimensions [Q, [7], ^ or view the Higgs 
boson as a Goldstone boson to explain its relatively low mass, i.e. little Higgs models [0]. A further 
possibility is a strongly interacting theory which might explain the mass of the W and Z-bosons 
even without a fundamental Higgs boson. There is of course another disturbing logical option: 
we see neither a hint of the Higgs boson nor any sign of BSM physics. Although this would be a 
nightmare for experimentalists it is actually an exciting option for theoreticians. Such a scenario 
can only be realised if something "invisible" disguises a low mass Higgs boson or if there are some 
strange interactions active in or beyond the electroweak sector which fake a light Higgs boson in 
the precision measurements by a yet unknown effect. Either we have to add an invisible sector to 
the Standard Model or we need to question our quantum field theoretical description thereof, both 
would be truly exciting ! 

The discovery potential of the Higgs boson has been investigated thoroughly in the experimen- 
tal studies of ATLAS and CMS [jK], 1 1 ] and signal significance plots like the one in Fig. [j] show 



'if it is ever possible to verify phenomenologically whether string theory is realised in nature is of course a different 
question. It is allowed to be sceptical here. 
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the importance and relevance of different search channels depending on the Higgs mass. Note that 
in most experimental studies K-factors, i.e. the effect of higher order corrections, have not been 
taken into account yet because higher order corrections to many background processes were either 
not known or the computer codes of known calculations are not flexible enough. Higher order cor- 
rections have been considered where available in recent experimental studies of the LHC physics 
potential [Q. 

Preferably discoveries at the LHC should not depend on the theoretical status of theoretical 
simulations. For example the diphoton decay of the Higgs boson should lead to a clear peak over 
the falling continuous background in the invariant mass distribution. This background should be 
measurable from the sidebands. Of course, a precise prediction is still needed for realistic studies 
before measurement [[1^, |l6|] and apart from that, who knows, whether the diphoton channel is not 
further contaminated by new physics? Other discovery modes, e.g. the H — ► W + VK~ discovery 
mode [17], has no clear signal and background separation and background subtraction has to rely 
on extrapolations of background to signal regions which do depend on the precision of theoretical 
input. 

In what concerns the evaluation of leading order estimates for LHC cross sections based 
on tree-level matrix elements the situation is quite satisfying as many tools have been devel- 
oped in the pre-LHC era, e.g. Alpgen @, COMPHEP [Q, Amegic++ @, COMIX [||], 
FeynArts/FormCalc Q, GRACE [|23|], HELAC [||, ||,||], Madgraph/Madevent [||, ||], 
Whizard pgh, and most of them are publicly available. However, any partonic fixed order cal- 



culation is of limited use for experimental studies which need realistic hadronic final states. This 
is the reason why merging of partonic amplitudes with partem showers has been an active field 
recently. The standard Monte-Carlo tools which contain partonic showering and include a hadro- 
nisation model are currently Pythia [30], Herwig [ ]3~ij ] and Sherpa [32]. Evidently, the more 
hard matrix elements are included (or better can be recursively evaluated) in these tools, the more 
accurate is the description of shapes of distributions and jet structure. In this respect Sherpa is 
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especially promising, as it is based on the matrix element generator Amegic++ [20]. On the other 
hand, the overall normalisation is predicted only very unreliably as long as virtual higher order 
corrections are not taken into account. This brings us to the main topic of the talk which is the 
computation of these higher order corrections. 



2. Framework for QCD computations 

The evaluation of production cross sections in hadronic collisions is based on the factorisation 
property of QCD. Whenever the hard scale of the process is significantly higher than the hadronisa- 
tion scale of around one GeV a differential cross section can be obtained by summing incoherently 
over the initial partons inside the hadron and can be decomposed into long and short distant contri- 
butions. 

dc(HiH 2 -> <pi + . • . + <Pn + X) = Yj] dx x dx2 fj/ Hl (xuHf) fi/H 2 (x2,Hf) 

j,i J 

x ( ia(parton J (xiPi)+parton / (x 2 P2) -> <pi + • • • + 9n, a s .(ju), ju f ) 

Here <pi + . . . + (p^ stands for some partonic Af -particle final state. The parton distribution functions, 
fj/ H (x,fip), which stand for the probability to find a certain parton in an incoming hadron, have 
to be determined experimentally at a certain scale. Their scale dependence on the other hand is 
determined perturbatively through the DGLAP-evolution equations. 

The factorisation scale {pLp ) dependence is stemming from the infrared (IR) structure of initial 
state singularities. In contrast, the renormalisation scale (Hr) dependence of the strong coupling 
constant is an ultra-violet (UV) effect. A hadronic cross section in leading order (LO) perturbation 
theory is thus plagued by a logarithmic dependence on these scales. This dependence is aggravated 
when the number of coloured partons increases due to higher powers of a s . A LO computation 
in QCD has thus only limited predictive power. By including next-to-leading order (NLO) correc- 
tions, the leading logarithmic dependence cancels and one obtains a much milder dependence on 
the unphysical scales. 

The evaluation of NLO corrections contains a loop and real emission part. Schematically 



^LO: 3BEr ^NLO, virtuaL ^)Qr ^NLO,reai: ^fnC 

The loop contribution has the same particle content and kinematics as the LO contribution. 
The real-emission contribution is tree-like but contains an additional parton. Apart from initial state 
collinear divergences, which are absorbed by the parton distribution functions, the IR divergences 
cancel between the virtual corrections and the integration of the real emission contributions over 
the soft/collinear phase space regions. In actual calculations the IR divergences are typically treated 
by subtracting terms which have the same IR limit as the real emission corrections and do allow 
for the analytic soft/collinear phase space integration. Adding back the integrated terms allows 
to represent the 2 — > N (NLO virtual) and 2 — > N + 1 (NLO real) part as independently finite 
contributions, which is mandatory in numerical implementations. 
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O = Olo + Gnlo 
Olo = J dPS N — \M 

Onlo = / dPS^ — a s 



LOI 



.c/t 



to<LO, virt. +Ko^O, virt. +1 / 



+ / dPS N+1 — a s 



I^NLO,reall 2 -L' 



(2.1) 



The most wildly used variant of these subtraction methods is the so-called dipole subtraction for- 



malism developed by Catani and Seymour [33] which has been generalised to the case of massive 
partons in p4|]. We see that the ingredients of a NLO QCD calculation consist of (i) tree level 
amplitudes, (ii) one-loop amplitudes and (iii) dipole subtractions terms, As was pointed out 
above, many automated tools are available for the evaluation of the tree-level contribution. More- 
over, different groups have recently implemented the dipole formalism in a process independent 
way [QUIZ, 



], and public codes exist which help to perform this labour intensive part of an 
NLO computation. The real bottleneck of an automated approach for NLO QCD computations is 
thus the evaluation of the one-loop matrix elements. This will be the main subject of the remainder 
of this talk. 



3. New methods to compute one-loop amplitudes 



3.1 Overview of achieved results 

The complexity of loop calculations grows rapidly with the number of external legs. What con- 
cerns computations with two initial and two final particles, basically all relevant LHC processes are 



evaluated meanwhile at NLO in QCD and public codes, like for example MCFM [f39|, |40j], the PHOX 
programs for observables including photons pT]], and VBF0NLO for weak boson fusion processes 
p2[], are available for experimental studies on the partonic level. Methods to combine fixed order 
calculations with partem showers have been worked out and implemented for many important pro- 
cesses, see for example MC0NLO [Q, POWHEG @, GRACE [||] and ©0]. Also the SANC 



3J, POWHEG pj, GRACE and [_ 
collaboration has promoted their framework to include QCD processes and progress has also been 
reported during this conference [^, ^9, 50]. A relevant gap of 2 — > 2 processes has been closed re- 
cently. The off-shell vector boson pair production cross sections for gg — > W* W* , , y*Z* , Z*Z* 
are available now as public codes GG2WW [^T], 52] and GG2ZZ [53]. These processes describe a 
relevant background component to H — > WW* /ZZ* also below threshold. These processes are in- 
duced by loop diagrams but are conceptually of LO type. Technically one has to square one-loop 
amplitudes which is numerically more delicate than computing interference terms between tree and 
loop contributions. 

Adding an additional final state particle increases the complexity of one-loop computations 
already considerably, and before 2005 only a limited number of 2 — > 3 processes had been evalu- 
ated and implemented in computer programs. NLO predictions have been existing for pp — ► jjj 
[p4|], pp —* yyj IP3> 01' PP ^ H jj in gluon fusion [|57|] and weak boson fusion (WBF) [p8l B^], 
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and pp —>■ Htt [J60, 61]. Between 2005 and this date most of the relevant 2^3 processes for LHC 
have been completed and the list of accomplished tasks contains also the processes pp — > HHH 
@i, pp W// in WBF 2 [g || ^], pp 

@ 01, -> VVV @ ^|] and /j/j -> ttZ, bbZ [0, ^|. For the important pp 
also electroweak 



ff/||6g, |69|], pp -> WWj 
Hjj process 



7|, [78p and sub-leading interference terms have been evaluated [ |79| , |8 
A discussion of the phenomenological importance of the different processes can be found in [p3]]. 
Note that the evaluation of the relevant amplitudes is only about half of the cake. The implemen- 
tation of a full NLO cross section into a reliable computer code is also a considerable effort due to 
numerical issues which will be discussed below. Further note that most of the corresponding codes 
are not public. 

The time-scale of all these calculations were of the order of person-years and it is no wonder 
that only a very limited number of results which deal with NLO 2^4 LHC processes can be 



found in the literature[84, 85 



], |87[ |88p. In the remainder of this section I will focus on the recent 
developments in this direction and report on progress which has been achieved by different groups 
and different methods recently. The two main approaches are based on either Feynman diagrams 
or on unitarity cuts. Before discussing these two approaches separately let us shortly comment on 
two important commonly used techniques. 

3.2 Colour and helicity management 

When evaluating gauge boson amplitudes or dealing with massless fermions it is very useful 
to work with spinor helicity methods, for an introduction see [f?9|]. It is well known that tree-level 
amplitudes are most efficiently represented in this way. The same is true for loop amplitudes. A 
massless Dirac spinor is already defined by two helicity states j^) defined by 



o 



n ± n(jfc) , (k ± \ = v(k)n T 



EL 



Helicity amplitudes can be written in terms of spinor products which are complex numbers: 

(kq) = (k-\q + ) , [kq] = (k + \q-). 
Massless gauge bosons like gluons and photons can be expressed using the same building blocks 



[kq] 



By construction one works with the two physical degrees of freedom. The auxiliary vector q 
defines an axial gauge. In full amplitudes its dependence drops out due to gauge invariance. The 
generalisation to massive particles is also well-understood. A public implementation of spinor 
helicity methods can be found in [j9(J. 

Once the helicity amplitudes are evaluated, the squared amplitude can be obtained numerically. 
If there are K helicity amplitudes labelled by {Xj} one only has to evaluate K products 



si 



Xj 



(3.1) 



"The six external particles are connected via at most 5-point one-loop topologies. 



6 



LHC phenomenology at NLO QCD 



otherwise K 2 terms have to be evaluated. To give an example, the TY-gluon amplitude can be 
decomposed into 2 N helicity amplitudes. Due to parity invariance only 2 N ~ l have to be evaluated. 
Nonetheless we see an exponential growth of terms in N. Another source for exponential growth 
in N-point amplitudes is the colour structure. A widely used colour decomposition is based on the 
following two elementary rules valid for an SU (Nc) algebra with generators T a in the fundamental 
representation: 

;rabcpc rparpb T>brra 

l J 1 ik — 1 ij 1 jk ~ 1 ij 1 jk 



T£Tjl = -(SiSl-—8i8/ 



Any colour structure can be mapped to simple Kronecker deltas which indicate the colour flow. 
This colour flow decomposition, described for example in [^l], ^], is going back to 'tHooft who 
applied a double line notation for gluons in the context of the 1 /Nc expansion. Any amplitude 
with Np fermion lines and Nc gluons can be decomposed in terms of AM = (Nf +Ng) ! products of 
colour flow lines. The different elements are labelled by the permutation group Sn acting on colour 
line indices. 

£ <K) , \c a ) = 8^8^...8^ (3.2) 

aeS N 

For gluon amplitudes this leads to an over-counting of independent colour states, as some singlet 
contributions are actually forbidden. The N-gluon amplitude has (N — 2) ! independent colour 
states. Still, asymptotically we notice a factorial growth of terms for large N. A way to fight 
the exponential growth both in helicity and colour components may be to resort to Monte Carlo 
sampling over these degrees of freedom, see for example [^] . 

3.3 The unitarity based approach 

The idea to use the analytic structure of scattering amplitudes to determine their explicit form 



is very old [ ]94| , |95Q . The well known non-linear relation between the transfer matrix, T, implied by 



the unitarity of the S -matrix, 5=1+ iT, 

S^S = 1 => 2 Jm(T) = T^T , (3.3) 

leads to a relation between the imaginary part of one-loop amplitudes and sewed tree-level dia- 
grams, schematically: 

Here the right hand side stands for the sum over all two-particle cuts. The latter correspond 
to squared tree-level amplitudes integrated over the respective two-particle phase space. Using 
the standard Feynman diagrammatic approach it can be shown that any one-loop amplitude with 
massless internal particles can be decomposed in terms of known scalar integrals with two to four 
external legs 3 4=2,3,4, and a rational part, ffl, which does not contribute to the imaginary part of the 




If masses are present one also has to consider one-point integrals. 
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Figure 2: Multiple cuts can be used to fix integral coefficients of amplitudes. 



amplitude. The coefficients of the integrals, C k , and the rational term are rational polynomials in 



terms of spinor products, see Eq. (3.1), and/or Mandelstam variables and masses. The imaginary 
parts of the different scalar integrals can be uniquely attributed to a given integral. Subsequently 
unitarity implies that the knowledge of the imaginary part of the amplitude defines the coefficients 
of the scalar integrals. 

^l-loop= £ C k I k +& Im^_ loop = £ C k \m{I k ) (3.4) 

£=2,3,4 fc=2,3,4 

Note that there are many different scalar two-, three- and four-point functions present in a given 
process. This is not indicated in this schematic equation. We see that tree amplitudes and phase 
space integrals in D=4 dimensions already fix a large part of the amplitude without ever being 
forced to evaluate any one-loop diagram. By turning the argument around we learn that one-loop 
diagrams are to a great extent determined by the pole part of the internal propagators which are 
nothing but cut propagators. The rational part is actually stemming from the UV behaviour of 
the one-loop amplitude, it comes from terms (D — 4)1® ~ 1. In amplitudes with an improved UV 
behaviour, such as susy amplitudes, these terms are absent, and this explains why these amplitudes 
are fully determined by the D=4 cut structure. This fact has seen many applications in the last years, 
see e.g. [fMj). As the rational part is not determined by the D=4 cut structure, one has to resort to 
more sophisticated methods, like for example on-shell recursion relations [ p7| ] or D-dimensional 
unitarity [^]. 

This method has been developed and applied to a great extent by Bern, Dixon, Dunbar and 



Kosower in the 90s [99, 100], for a recent review see [96]. It has seen a renaissance after the Santa 



Barbara workshop in 2003 where Witten pointed out a relation between certain string theories and 



QCD amplitudes using twistor methods [101]. The revived interest lead to new efficient recursive 



evaluation procedures for multi-parton amplitudes [102, 103, 104]. A main feature of recent vari 



ants of the method is that multiple cuts are used to determine integral coefficients, which goes under 



the name of generalised unitarity [102, 103]. Cutting four lines in an N-point topology amounts 
to putting the corresponding four propagators on-shell. This procedure fixes the associated loop 
momentum completely and the coefficient of the related box diagram, see Fig. ||(a), is given as a 
product of tree diagrams. 

Ro BW,v =±> r - 1 V aM 4 ( 3 ) 
ng-l^fa;. =?• c 4 — A tree "nree 71 tree A tree 

z a=± 
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Figure 3: The histgrams of the evaluated events are plotted over the logarithmic relative error 
log ( | tumeric — MargeOI/IMargetl) °f the pole and finite parts of the MHV N=6,7,8-amplitudes. The 
shaded region (blue) relates to the finite contribution, The unshaded full line region (red) to the l/e pole part 
and the dashed to the double-pole part [113]. 



One has to sum over the two complex solutions of quadratic on-shell conditions. This fixes the 
loop momentum. An important point is that the tree level amplitudes already incorporate gauge 
invariance manifestly. In Feynman diagram computations many graphs have to be combined to 
result in a gauge invariant expression. The power of this method has been demonstrated by the 
analytic evaluation of the cut-constructable part of the six-gluon amplitude [104]. Combined efforts 
of many groups were needed to compute the different pieces of this amplitude, see [105, 106] for 
a collection of all formulae. The rational part has been first provided by a Feynman diagrammatic 



computation [ ]107| , |10q , |109j]. Note that the six-gluon one-loop amplitude is part of a full NLO 
computation of four jet production at the LHC and thus of phenomenological interest. There exists 
also a very compact result of the six-photon amplitude which was obtained using cutting methods 
[11C, 111 ]. As has been shown in [1T2], the rational part is zero in this case. 

The unitarity method, using multiple cuts, has recently been implemented in a numerical code 
by the Blackhat collaboration [ ]113| , |114 ] and by Ellis et. al. [115] (discussed below). To 
this moment the Blackhat collaboration has provided two remarkable applications, the colour 
ordered part of the (N=6,7,8)-gluon amplitudes in a special helicity configuration, H--+ anc j 
the leading colour contribution of qq — > Vggg (this includes crossing related amplitudes) [ ]1 16| , pfj \. 



Important information on the numerical behaviour of the method could be gained in [116]. 
The large number of algebraic operations in these complex amplitudes typically lead to round 
up errors in certain exceptional phase space regions. Any automated method has to come up 
with a reliable fail-safe procedure if the numerical cancellations are spoilt by the finite number 
of computed digits. In Fig. ||]a comparison of the MHV 6,7,8-gluon amplitudes between the nu- 
merical implementation and the known analytical result is shown. The logarithmic relative error, 
log (I tumeric — target) I /I target I)' i s plotted for the double, single-pole and the finite part of 
these amplitudes for 100.000 phase space points in Fig. |3| While the bulk of all points is evaluated 
reliably using standard double precision, it can be observed that the distribution has a tail which for 
a few percent of phase space points indicates a precision loss of almost all digits. Such a situation 
typically leads to an unstable behaviour of adaptive Monte Carlo event generators, which tend to 
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Figure 4: The leading colour NLO prediction for pp — ► V jjj compared to Tevatron data. The leading colour 
prediction is expected to match the full result within a few percent. As expected the error bands due to scale 
variations reduce drastically when moving from LO (brown) to NLO (grey) [p7n. 



sample points exactly in regions which induce a large variance of the result. In the given plot one 
observes that the tail is cut at a precision of around 10~ 3 to 10~ 4 . Beyond that point multi-precision 



libraries Ql 17Q are used to avoid a dangerous loss of precision. 

The evaluation of the loop amplitude to qq' — ► Wggg is the necessary ingredient for a pre- 
diction of the important Standard Model background pp — ► Vjjj. Such multi-jet plus lepton plus 
missing energy signals occur for example in supersymmetric extensions of the Standard Model. 



The Blackhat collaboration has confronted their leading colour prediction for pp — ► Vjjj [ |87[ ] 
with Tevatron data, see Fig. ||. As expected the stability of the prediction under scale variations 
is largely improved. The leading colour approximation is expected to be phenomenologically jus- 
tified within a few percent for Tevatron kinematics. Technically this approximation simplifies the 
calculation considerably, as many topologies are suppressed and can be neglected in this way. 

In 2006 Ossola, Papadopoulos and Pittau (OPP) have proposed another unitarity approach 
which was based on the question how to reconstruct one-loop amplitudes from the integrand of 



Feynman integrals [118]. Let us consider a sub-amplitude where all external particles are ordered. 
The full amplitude can be obtained by adding up different permutations. The leg-ordered amplitude 
can be obtained after summing the corresponding Feynman diagrams written such that the sum is 
under the integral and has a common denominator of propagators D\ . ..D^ with Dj = (k — rj) 2 — 

2 



<\ V \ 



dPk- 



D { ...D N 

The fact that general one-loop amplitudes can be written in terms of at most 4-point scalar 
integrals implies that the numerator can be written schematically in the following way 

JT(k)~ £[C 4 + C 4 (*)] f[Dj+ £ [C 3 +C 3 (k)} fl D J 

boxes j^box triangles j^triangles 

N 

+ £ [C 2 + C 2 (k)] n D J + --- < 3 - 5 ) 

bubbles j^bubbles 
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In this expression the coefficients C4, C3, C2, . . ., correspond to the actual amplitude coefficients of 



the various box, triangle and bubble integrals defined above, see Eq. ( |3.4| ). While these coefficients 
are independent of the loop momentum, the numerator function also contains spurious k-dependent 
pieces, C4, C3, C2, . . ., which have to vanish upon integration over the loop momentum k. OPP asked 
the question how these spurious coefficients can be defined without doing the actual integration. 
They concluded that one simply has to evaluate the on-shell conditions of the various propagators 
which is in one-to-one relation to the unitarity cut methods using multiple cuts. This approach 
defines the coefficients of the terms C4 ,€3, C2, ■ ■ ., as polynomials in k. The coefficients of the box- 
, triangle-, bubble terms are defined by respectively 2,7,9 /^-polynomial coefficients. OPP proposed 



to extract these coefficients by numerical interpolation of the various polynomials in k [1 18]. They 



also define a procedure to access the rational part of the amplitudes [ 1 19, 120Q. This method works 



also for individual one-loop Feynman diagrams and it has been implemented in a public computer 



code, CutTools [121]. The method has been applied successfully to the evaluation of six-photon 



amplitudes with massless and massive internal fermions Qlll[ |. The method has also been applied 
to the QCD corrections of triple vector boson production at the LHC. Here the amplitude has been 
integrated to obtain the total and differential cross sections for relevant observables [[72|] . 



A variant of these numerical unitarity based methods has also been investigated in [115] and 



has been extended to a D-dimensional approach [ |122[ , |123| ]. In the latter case also 5-point functions 
have to be included in the function basis. The key observation is that by adding higher dimensional 
information the rational part of the amplitude is also determined. A D-dimensional numerator 
function of a leg ordered amplitude is of the form 



jVd = ,yV^ + (D — 4)^_ 4 = £ [C 5 + C 5 (k)} ]J Dj + [r.h.s. of Eq. Q] 

pentagons pentagon 

This approach guarantees that the evaluation of the rational part is now on the same footing as the 
cut constructable part. The method relies on the evaluation of amplitudes in dimensions different 
from four. It can be viewed as a D-dimensional generalisation of the OPP approach. The authors 
observed that the knowledge of tree amplitudes evaluated for complex momenta which are defined 
through the on-shell conditions of propagators, or equivalently through multiple cuts, defines the 
complexity of the algorithm modulo the number of different possibilities to cut a given N-point 
topology. The latter scales like ~ N 5 , as any set of five propagators have to be put on-shell in a 
given diagram. The recursive evaluation of tree-amplitudes scales approximately like ~ N 4 . Thus 
it can be expected that a colour ordered multi-parton amplitude can be evaluated by such a polyno- 
mial complexity algorithm, or short ^-algorithm, [ 122| , 123 ]. As was pointed out above, for full 



amplitudes one needs to sum over different helicities and colour structures which unavoidably turns 
on an exponential growth and thus one ends up with an exponential or $ -algorithm for full ampli- 
tudes. Again Monte-Carlo sampling over colour and helicities might help to single out the relevant 
contributions numerically in a more efficient way. If a polynomial algorithm can be realised in this 
way remains to be shown. 

The potential of this valiant of the unitarity method has been demonstrated by evaluating single 
phase space points of different helicity amplitudes of the colour ordered Af-gluon amplitude for up 
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Figure 5: Evaluation time of colour ordered A^-gluon helicity amplitudes. The dependence of tree- and 
one-loop amplitudes on N confirms the expectation of a polynomial growth proportional to ~ A^ 4 and A^ 9 



respectively [124] 



to N = 20 legs | |124[ ]. Sampling over a large number of phase space points shows a good numerical 
performance concerning speed. 

In Fig. H the polynomial behaviour of tree and one-loop amplitudes is shown as a function of 
N. The polynomial behaviour of the tree, ~ N 4 , and the one-loop amplitudes, ~ N 9 , is confirmed 
for the multi-gluon amplitudes. At the moment it seems not feasible to produce such a plot based 
on Feynman diagrams ever, as gauge cancellations and colour decomposition are not manifest 
and one needs the evaluation of the order of AM form factors. An independent implementation of 



the same method and amplitude, confirming the result, has been presented recently fll25| , |12q ]. 
Other applications of the method are the evaluation of the amplitudes gggtt — ► [ ]123| ] and the 
leading colour contributions to qqWggg, qqq'q'Wg J86|, 127]. The latter are also implemented in 



the Rocket code [ ]124| ] which handles numerically problematic phase space points by using higher 
precision libraries. Again, ultimately multi-precision libraries are in use to guarantee numerical 



reliability [ ]117| ]. 

From the fast progress made in the last two years it can be expected that further cross section 
calculations for LHC processes with four or more partons/particles in the final state will become 
available using the discussed unitarity based methods. Note that apart from the multiple cut meth- 



ods described so far, also single cut techniques are under development [ |128Q . 



3.4 The Feynman diagrammatic approach 

The traditional method of performing amplitude calculations is by starting from Feynman 
diagrams. It is based on representing a one-loop amplitude as a sum of diagrams sorted by colour 
structures 

■</' A il>l-'"P = I f {Cl} ^ } 
{a},a 

-J iK D/2 D\...Dx -L^u-.H* h (Pj^j) 
d D k k^...k^ 



n \Ph m i) = J iK D/2 Dl ^j) N > D j = ( k - r j) -mj,rj = Pl ...+pj. 



(3.6) 
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Figure 6: Left: Comparison of the scale dependence for the LO (dashed) and NLO (full) cross section 
pp(qq) — » bbtt at LHC. Right: Effect of an invariant mass cut on the b-pair. The scale variation is indicated 
by the two bands, LO (blue), NLO (red) [fl. 



The kinematical functions can be represented by tensor one-loop integrals. The latter can be evalu- 



ated recursively in momentum space using the Passarino-Veltman method Q129j , |130p . Alternatively 
one can map the momentum integrals to Feynman parameter integrals which also can be reduced 
recursively [|3l], |33|, [134], |l3|, |L37] ]. 



^■■^(r jl ,...,r jr ^I^ 2m Uu--.Jr) 

f°° ,\j „ , . ^-i \ Zh---Zi 



N 

L 

1=1 



{-\z-y-z) N - D ' 2 



3>ii 



(r,- — rA —m^ — rn 2 



J ' 



The end-point of the different recursion algorithms are scalar integrals with no numerator structure. 



Public program libraries exist to evaluate the latter Q138| , |139| , |140| , [141] , |142| ], 

Over the years quite a few groups have gathered a lot of experience in efficiently implementing 
Feynman diagram calculations. Cross section calculations, with up to five external particles, can 
generally be mastered now, as is demonstrated by the long list of accomplished tasks discussed 
above. For 2^4 processes only a view NLO computations have been accomplished using Feyn- 
man diagrams, e + e~ 



/1/1/2/2 [MJ, e + e- HHw [|144J, and 77 bbtt [Jl45|]. Processes of 
this complexity relevant for LHC phenomenology are under construction right now. Progress on 
the computation of the important LHC process pp — > bbtt was reported in 2008 [84], this group 
presented the full cross section computation of the quark induced subprocess qq — > bbtt at NLO in 
a s . The amplitude can be written in terms of rank three 6-point tensor integrals 



sf(qq -> bbtt) 



E r jhhh 
1 /:/''\^ 6 



(3.7) 



In this calculation the algorithm described in [ |146[ ] was used for the tensor reduction of 5- and 
6-point integrals, otherwise Passarino-Veltman reduction was used. Note that the gluon induced 
subprocess gg — ► bbtt has a more complicated tensor structure, nonetheless a full result for this 
important Standard Model background can be expected in the near future. 

Another collaboration is currently pressing for an automated, Feynman diagrammatic evalu- 
ation of multi-leg one-loop amplitudes. The GOLEM collaboration, where the acronym stands for 



13 



LHC phenomenology at NLO QCD 



General One-Loop Evaluator of Matrix elements, set up an automated reduction framework which 
takes special care of numerical issues by providing alternatives for the evaluation of one-loop in- 



tegrals in critical phase space regions [ |135| , |147| , |141| ]. In this approach rank R N-point functions 
are reduced first to 6-point tensor integrals and the latter are then expressed by Feynman parameter 
integrals like 



jD,D+2, ■ 
*N=3A Ul) 



1 4 



n<M(i-i>) 



(=1 



7=1 



(-h-y-z-idy- D / 2 



Spurious numerical problems due to so called inverse Gram determinants are avoided by having 
the option to evaluate the form factors in two different ways, either by numerical means or by 
reducing the form factors numerically to a scalar integral basis. The method is designed to allow 
for an efficient isolation of IR and UV divergences from the result and as such very well adapted 
for QCD calculations. The method can be also used to evaluate the rational part of an amplitude 



only [ 112Q . All 6-point form factors are coded in a f ortran 95 code golem95, which is the first 



public library for such form factors Q141[ ]. The code golem95 relies on the evaluation of scalar 
integrals. In the present version only integrals with zero internal masses are implemented. Another 



implementation of 6-point form factors has been presented recently Jl48| , |149| ]. 

A full amplitude evaluation relies on diagrammatic input. The GOLEM collaboration uses pub- 



lic tools like QGRAF [ ]150| ] and FeynArts [ |22[ ] for this step. The colour algebra and the helicity 
management is done as outlined above. Subsequently two independent strategies are followed. 
Firstly, the diagrammatic input is converted to a form factor representation 



■cK +2 » D+4 (u^)+c^» D+2 (u^) + . 



(3.8) 



exported to a f ortran code and then linked to the form factor library golem95. Secondly a 



fully symbolic representation of the scalar integral coefficients is generated using FORM [151] and 
MAPLE 



(3.9) 



Using these implementations several computations of relevance for the LHC have been performed 
[152, 52, ^3[ |8l| ]. The viability of the approach for 2 — > 4 processes was shown by the evaluation 
of the 6-photon amplitude mediated by a massless electron loop [ |1 Uj ). This amplitude has been 
used as testing ground for various methods and perfect agreement has been obtained in all cases 
JT53|, [TTT1, p4] , [T55l ]. 

Currently the NLO corrections for the process pp — » bbbb are under construction using the de- 
veloped computational tools. Because the 4b-quark background needs to be known as precisely as 
possible in the context of two Higgs doublet models, this process was included in the Les Houches 
"experimentalist wish-list" [S3]. The amplitude consists of two different initial states qq — > bbbb 
and gg — ► bbbb. The first can be represented by about 250 Feynman diagrams among which 25 
pentagon and 8 hexagon diagrams can be found, see Fig. [7] for a selection of LO and NLO dia- 
grams. Note that only the pentagon and hexagon diagrams are computationally challenging. The 
other topologies are relatively simple and do not present a problem apart from book keeping which 
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Figure 8: Comparison of double (full) vs. quadruple precsion (dashed) evaluation of the local K-factor 



function and the single/double pole zero coefficient as defined in the text [88 1 



is not an issue in automated approaches. The total number of diagrams is thus not a good mea- 
sure of the complexity of a calculation. All one-loop Feynman diagrams have been evaluated in 
two completely independent ways, by symbolic reduction to scalar integrals and by using the form 
factor decomposition in combination with the golem9 5 library as discussed above. The virtual 
&((X S ) correction is a interference term between the LO tree and NLO loop amplitude. After UV 
renormalisation, IR divergences in form of 1/e single and double poles remain. The IR structure 
of one-loop amplitudes is well-known and a finite function can be obtained after adding an ad- 
equate IR compensation term, as defined for example by the Catani-Seymour insertion operator 
(ja4o|I(e)|^Lo) [Hi- In this way one can define a finite local K-factor function 



K 



/ LO | 2 + 2Re(^ 



WLOMrt ) 



UV/IR subtractions 



(3.10) 



To investigate the numerical performance of the approach, 200.000 random phase space points 
have been evaluated using double and quadruple precision, see Fig. ||. In the figure the numeri- 
cal precision of the evaluation of the local K-factor and the cancellation of the IR single/double 
pole are compared. The result indicates that the size of the finite K-factor is a good indication of 
the numerical accuracy of the evaluation. Large K-factors are mostly induced through numerical 
round-up errors. The cancellation of single and double poles, which is used by other collaborations, 
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Figure 9: Renormalisation scale dependence of the LO cross section (full red) and including the UV/IR 
subtracted virtual correction term (dashed black) [|88||. 



can also be applied to decide about the quality of the evaluation but seems to be less indicative. In 
Fig. |9] the improvement of the scale dependence is shown, if the subtracted virtual correction term, 



as defined in Eq. (3. 10), is included. The cross sections have been evaluated using the experimental 
cuts T] < \2.5\, AR > 0.4, p T > 25 GeV and the NLO pdf set CTEQ6.5. A further improvement 
of the behaviour is to be expected, because the logarithmic factorisation scale dependence is not 
yet compensated. To obtain the full NLO result one has to add the dipole subtracted NLO real 
emission corrections. For the full process, pp — > bbbb, the gluon induced sub-process also has to 
be included. These parts of the calculation are in progress. 

Note that the evaluation strategy of the virtual NLO corrections differs from the conventional 
approach of integrating a NLO cross section directly. In the given example the LO cross section 
has been evaluated first to produce an unweighted event sample. The subtracted virtual correction 
term defines a finite K-factor function which is simply used to reweight these LO events. In this 
way interference effects between adaptive Monte-Carlo integration and numerical round-up errors 



are completely avoided. For a more detailed discussion see Q85Q. A similar strategy has also been 



applied in Q67, 72]. 



The goal of the GOLEM collaboration is to provide an automated framework for one-loop am- 



plitude computations. The GOLEM implementation is sketched in Fig. 10. The input files which 
contain process and model dependent information are processed with a python script which con- 
trols QGRAF and FORM codes. The diagrammatic output, including UV and IR subtraction, is 
defined in terms of form factors and written to fort ran 95 code. The latter is linked to the form 
factor library golem95. At the moment the code can be compiled in double or quadruple preci- 
sion but in the future a more dynamical solution will be adopted. The resulting computer program 
provides a function to reweight LO event samples stemming from tree level event generators. For 



the shown results Whizard [29] has been used 



3.5 Is there a preferred method? 

We have seen that both the unitarity based and also the Feynman diagrammatic method are 
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Figure 10: Flow chart of the GOLEM implementation [B8H. 



sufficiently well developed meanwhile to allow for cross section calculations with more than five 
external legs. While the unitarity based methods were initially used to obtain compact analytical 
expressions for helicity amplitudes, the quest for automation lead recently to various numerical 
implementations. Higher precision libraries play a prominent role in the related computer programs 
which shows that any given method seems to be doomed to deal with exceptional phase space 
regions in one way or the other. Higher precision libraries seem to be a natural cure for numerical 
problems. As the number of problematic phase space points is typically only a few percent of all 
points, increased evaluation times are not a severe issue 4 . 

Different methods will soon provide cross section predictions for partonic 2^4 processes rel- 
evant for LHC phenomenology. As was explained above the scaling properties of one-loop N-point 
amplitude computations seem to indicate that the unitarity based methods are clearly preferable. 
The factorial growth in diagrammatic computations has to be compared to polynomial algorithms 
now. For example, the evaluation of a colour ordered multi-gluon helicity amplitudes scales like N 9 
whereas a Feynman diagrammatic calculation would involve about 2 N leg ordered diagrams with 
roughly T(N) = (N — 1)! form factor evaluations each. However, let us ask the question, whether 
the asymptotic behaviour gives a guideline which method is preferable for phenomenological appli- 
cations. The logarithmic ratio of the asymptotic behaviour, \og(N 9 / (T(N)2 N )) is plotted in Fig. 11. 
The figure indicates that for phenomenological relevant multiplicities, say N <8, the asymptotic 
scaling behaviour is not a good measure for the question which method is preferable. Furthermore 
the true complexity of a method should be measured with full processes, taking into account all he- 
licity amplitudes, relevant colour structures and sub-processes. It is important to note in this respect 
that especially for background cross sections severe experimental cuts will be applied. Experience 
shows that any approximation or assumption might be invalidated in that way. For the full task no 

4 As discussed during this conference a hardware implementation of quadruple precision in processors is not yet in 
sight although technical standards have been already defined for such a step. The particle physics community would 
surely be grateful for such a development. Unfortunately fundamental science issues do not play the role they deserve in 
business plans of computer chip producing companies, although most technological developments are essentially based 
on progress in fundamental science! 
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Figure 11: The logarithmic ratio of the asymptotic behaviour, \og(N 9 / (T(N) 2 N )), is plotted versus N. 



^-algorithm exists. 

For ,/V < 8, implementation issues, like efficient caching and the re-use of recursive information 
can still be improved in Feynman diagram calculations. Together with the increasing computer 
power 5 it can be expected that both methods will deliver what they promise, namely cross section 
predictions for LHC phenomenology at the next-to-leading order level. 

At the end of this section I would like to add the comment that fully numerical approaches 



for Feynman diagram calculations have also been investigated by several groups [ |156| , |157| , |158 



153| , 159 , 155 ]. Given the huge world-wide computing resources any reliable integration method of 
multi-parameter integrals containing threshold singularities could lead to a solution of the one-loop 
problem without using any algebraic overhead. At the moment these methods are not competitive 
with the discussed approaches but there is certainly room for large improvements. 



4. Conclusion and Outlook 



In this talk I have argued that the LHC needs and deserves an effort to predict prominent sig- 
nal and background processes at the next-to-leading order level in QCD. Especially multi-particle 
amplitudes have a high power of a s (^l) which induce large scale variations. Absolute rates thus 
can not be predicted reliably with leading order Monte-Carlo tools and eventually this will ham- 
per the understanding of LHC data and the discovery of new physics. Many relevant Standard 
Model processes with three final state particles are meanwhile available in the literature beyond 
the leading order. Unfortunately the results are not always available as a public code which would 
be most beneficial for the experimental community. The theoretical community is at the moment 
enormously active to go beyond the given level of complexity and NLO predictions for processes 
like pp — > jjjj, pp —>W jjj, pp — > WW jj, pp — > bbtt, pp — > jjtt, pp — > bbbb, which are highly 
relevant for various Higgs and BSM search channels, are under construction. Many different meth- 
ods have been developed and are applied for the virtual corrections. Unitarity based methods look 

5 We learnt at this conference that the "multi-core" computer era will change to the "many-core" era were the number 
of idle and available processors will increase enormously. 
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very promising in this respect but for the required level of complexity both, the Feynman diagram 
and the unitarity based approaches, will provide valuable results for phenomenology. The avenue 
to fully numerical approaches to one-loop amplitudes is still not fully explored, but given the ever 
growing computer power this might be a promising direction avoiding any kind of algebraic ma- 
nipulations of initial expressions. 

Most groups move at the moment towards automated approaches which will provide platforms 
to do many computations with the same set-up. To cure the problem of numerical instabilities 
multi-precision libraries seem to be the accepted way out now. It would be highly appreciable, 
if the community could agree on a standardisation of in- and output formats for NLO codes and 
would make computer programs publicly available. The structure of one-loop computations is 
indeed very modular and a minimal agreement on colour and helicity management and on passing 
IR subtraction terms, which are basically universal anyhow, could result in transportable modules 
for virtual corrections. This would allow to use computations of different groups interchangeably 
by combining them with public tree level matrix element generators which should of course also 
contain the functionality of IR subtractions. In a next step the combination with partem showers 
could be obtained. Here the inclusion of NLO precision is also well understood, but again a focus 
on modularity and transportability would help to use synergies between different groups. 

The main conclusion of this talk is that the developments of the last few years are spectacular 
and point towards Monte Carlo tools at full next-to-leading order QCD level, at least as long the 
number of final state particles is not excessively high. The dominant and most relevant processes 
for Higgs and new physisc searches will certainly be available in form of flexible and reliable public 
computer programs and the analysis of LHC data beyond the leading order in a s will eventually 
become the standard. 
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